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Searching for zeroes: unconventional superconductors in a magnetic field 
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Abstract 

CJ 

> 1 | We review the results of the microscopic approach to the calculation of the anisotropy in the specific heat in unconventional 
O S uperconductors under rotated field. Treating vortex scattering on equal footing with the energy shift we find that the electronic 
^ specific heat may have minima or maxima when the field is aligned with the nodes, depending on the temperature and field range. 
^ We discuss the influence of the paramagnetic limiting and Fermi surface shape on the location of the inversion line. 
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1. Introduction 

One of the practical steps that help determine the ori- 
gin of superconducting pairing in novel materials is the 
determination of the shape of the energy gap. There are 
few techniques that directly probe the location of the ze- 
roes (nodes) of the superconducting gap in the bulk of un- 
conventional superconductors. Prominent among them are 
measurements of the specific heat and the thermal conduc- 
tivity under a rotated magnetic field. Here we review recent 
progress on theoretical underpinning of these experimental 
techniques, and their implication for realistic materials. 

The basic idea underlying the measurements is based on 
the understanding that in superconductors with zeroes of 
the gap on the Fermi surface (FS) the unpaired quasipar- 
ticles have momenta close to those of the nodal directions. 
Specific heat (thermal conductivity) measure entropy (its 
transport), and are therefore insensitive to the condensate 
of Cooper pairs. Consequently, if an experiment selectively 
probes excitations over a region of the FS, it is able to 
help locate the gap nodes. Spectroscopic techniques (e. g. 
ARPES) can resolve the gap structure in quasi-2D systems 
with large gap; however, in materials with low T c an alter- 
native method is needed. 

Magnetic field provides such an alternative. The field 
couples to the phase of the superconducting order param- 
eter, and drives supercurrents that, in turn, affect near- 
nodal quasiparticles. In the mixed state, when the field is 
only weakly screened, the interplay between the supercur- 
rents and single particle excitations occurs throughout the 



bulk, and is therefore accessible experimentally. 

In the semiclassical approach, valid at low fields, 
H -C H C 2, and low temperatures, T <C T c , the position- 
dependent superflow around the vortices, v s (r) shifts the 
energy of the unpaired electrons relative to the moving 
condensate, generating a finite field-dependent density of 
the quasiparticles. In the absence of the field the cost of 
a single electron excitation with momentum k is E^ = 
\/Ck + I Ak | 2 , where (k is the band energy, and Ak is the 
superconducting gap. In the presence of the Doppler shift 
due to slowly varying superflow this energy becomes E^ = 
.Ek — v s( r ) • k, and therefore some previously unoccupied 
states drop below the chemical potential and become pop- 
ulated. Since the Doppler shift outside of vortex cores is 
smaller than the maximal gap (supercurrent is below the 
pairbreaking value) , the populated states are located near 
the gap nodes in the momentum space. 

Pioneered by Volovik to predict that the density of 
states in a clean superconductors with lines of nodes varies 
as VH [T], this approach was extensively used to analyse 
nodal superconductors [2] • Early analysis aimed at high-T c 
cuprates, and started with the assumption of a cylindri- 
cal Fermi surface with vertical line nodes. Vekhter et al. 
considered the field rotated with respect to nodes in the 
conducting planes, and showed that the density of states 
(DOS) varies with the angle of field rotation [3]. Since the 
superfluid velocity is always in the plane normal to the 
applied field, the Doppler shift for the nodal quasiparti- 
cles, v s (r) • k„, depends on the angle between the field and 
the nodal direction, and vanishes when the field is along 
the node. Consequently, the DOS was predicted to have a 



minimum when the field is aligned with nodes, suggesting 
the measurements of the electronic specific heat as a test 
for the nodal directions. This behavior was indeed found 
in the vortex state of several unconventional superconduc- 
tors [4]. The results were qualitatively consistent with the 
predictions of the semiclassical theory, even though the 
experiments were done at higher T/T c > 0.1 and higher 
H/H c 2 than the corresponding range in the cuprates. 

Experiments on the anisotropy of thermal conductivity, 
k(H, T), under a rotated field continued for over 10 years, 
and unambiguously showed the superposition of the varia- 
tion due to difference in transport normal and parallel to 
the vortices ("twofold"), and the additional features due 
to nodal structure [5] . Theoretical analysis of these mea- 
surements remained incomplete until recently. The Doppler 
shift approach leaves the lifetime in the absence of impu- 
rities infinite, and therefore does not account for the scat- 
tering on the vortices. The magnetic field does influence 
lifetime via the self-consistency between the DOS and the 
impurity scattering, so that the field dependence of n is 
determined by the competition between the enhanced den- 
sity of states and the change in lifetime. Consequently, it 
is not clear whether rotating the field through the nodes, 
for a fixed direction of the heat current, gives minima (as 
the DOS) or maxima (due to reduced scattering). More- 
over, the connection between the local re(r) and the mea- 
sured value /c(H) for a given field direction is not fully es- 
tablished, and 'ad hoc' averaging procedures often entirely 
miss the twofold anisotropy. 

The need for a complete theoretical analysis, accounting 
for the superfiow beyond semiclassical treatment and ad- 
dressing the scattering on the vortices, was emphasized by 
comparison of the data on the specific heat and the thermal 
conductivity anisotropy in heavy fermion CeCoIns , where 
the opposite gap structure was inferred from the specific 
heat [6] and thermal transport [7j measurements. Such a 
description was recently developed by us [8|9j , and here we 
summarize the salient features of our approach, review the 
essential physics, and focus on a new aspect of the analysis, 
the effect of the Fermi surface shape. 

2. Microscopic approach 

Our approach is based on the quasiclassical formulation 
of the microscopic theory. In the spin and particle-hole 
space the matrix propagator, which depends on the direc- 
tion at the Fermi surface (FS), p , and the center of mass 
coordinate, R, and the mean field order parameter are 

\itr 2 l -g J \icj 2 A* J 

The equilibrium propagators obey 

[(e + ~Vir(p)A(R)) ? 3 - A(R, p) - a imp (TL; e),g(R, p; e)] 
+zv F (p)- V R g(R,p;e) = 0, (2) 



with the normalization condition g 2 (R, p;e) = — ir 2 l. 
In Eq.(2) vj?(p) is the Fermi velocity for the direction 
p, the vector potential A(R describes the applied field 
that we take to be uniform, and the impurity self-energy, 
tTi m p(R;e), is determined within the self-consistent t- 
matrix approximation. The order parameter is computed 
self-consistently with the separable pairing intercation 
y(P)P') = Voy(p)y(p'), where y(p) is the normalized 
basis function for the angular momentum representation. 
We consider quasi-two dimensional (open along z) Fermi 
surfaces and order parameters with vertical lines of nodes. 

We model the vortex state by the extended Abrikosov 
solution consisting of superposition of both the lowest and 
higher Landau level functions, <f> n (x), required for uncon- 
ventional superconductors, A(R) = J2 n A n (R\n) with 

R n = V CT ) —==■$„ )_ J v . 3 

For the field applied at an angle Oh to the z-axis, the 
anisotropy factor Sf — \J cos 2 Oh + («of| / v o±) 2 sin 2 Oh , 
where v 2 0X = 2{y 2 {p)vl l { Pz )) FS ,v 2 ^ = 2(y 2 (p)v 2 (p z )) FS , 
(•) FS denotes the Fermi surface average, i?n is the pro- 
jection of Vf on the z axis, and v±i with i = x, y is its 
projection on the axes normal to z. 

We use the Brantd-Pesch-Tewordt (BPT) approximation 
in replacing the normal (diagonal) part of the propagator 
with its spatial average, solve for the off-diagonal compo- 
nents in terms of g, and use the normalization condition 
to determine the quasiclassical Green's function. We then 
compute the density of states directly, and use the linear 
response theory in Keldysh formulation to determine the 
thermal conductivity. The justification of the approxima- 
tion, and the technical details of the approach are given in 
Ref. [9] for a rotationally symmetric Fermi surface in the 
shape of a corrugated cylinder. Here we briefly review the 
main results of that work, discuss its relation to the data 
on CeCoIns, and focus on the influence of the Zeeman field 
and the shape of the Fermi surface on the conclusions and 
the phase diagram. 

3. Results. 

In Refs. [819] we considered a Fermi surface given by p 2 — 
Px ^Py ~ i 7 " 2 Pf) cos (2sp z /r 2 p/), with s = r = 0.5, and the 
G?-wave gaps, y = v / 2cos2</), where <f> is the azimuthal an- 
gle at the Fermi surface. The unexpected result of Ref.[S] 
was the inversion of the anisotropic profile of the specific 
heat across much of the phase diagram in the T-H plane. 
In the low temperature, low field region the results of the 
microscopic calculation agreed with that of the semiclassi- 
cal approach, and the specific heat for the magnetic field 
applied along a nodal direction was lower than that for the 
field along a gap maximum. In contrast, upon increasing ei- 
ther the field or the temperature, the pattern changed and 
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Fig. 1. Origin of the inversion of the specific heat anisotropy. a) 
Low-energy density of states for H/H c 2 = 0.18 and field in the nodal 
and antinodal directions; b) Typical angle- resolved DOS (red) in 
zero field; c) Corresponding angle resolved DOS (red) under applied 
field. The order parameter profile is shown for comparison (green). 

the maximum, rather than a minimum of the specific heat 
was found when the field is aligned with the gap node. For 
a corrugated cylindrical Fermi surface reminiscent of that 
of CeCoIns, the inversion occurred already at T/T c ~ 0.2 
at low fields, and H/H C 2 ~ 0.5 at low temperature. There- 
fore over a very large part of the T-H phase diagram, away 
from the semiclassical low-energy regime, the theory pre- 
dicts the anisotropy of the specific heat opposite to that 
expected at low fields. As a result we suggested that the 
measurements in CeCoIns, performed at moderate H and 
T, were consistent with the d x 2_ y 2 symmetry. Our analy- 
sis of the thermal conductivity accounted for the first time 
both for the twofold anisotropy due to vortex scattering and 
for the structure due to nodes, and supported the d x 2_ y 2 
symmetry of the order parameter. 

The inversion is traced to the effect of vortex scattering 
as shown in Fig. 1. At low T the specific heat is related to 
the DOS N(p,e)/N = -Ini5(p,e) via 



C(T,H) = 



1 



+ 30 



f/t 



e 2 N(T,H:e) 
T 2 cosh 2 (e/2T) 



(4) 



and is most sensitive to the DOS at e ~ 2T. Inversion 
of the C(T, H) anisotropy at moderate T for a given H 
implies that the DOS anisotropy for the nodal and antin- 
odal field directions changes sign, see Fig. la): at low 
(moderate) energies N(H\\node,e) < N(H]\ antinode, e) 
(N(H\\node, e) > N(H\\ antinode, e)). Consider the angle- 
resolved density of states and recall that in the absence 
of the field the major contribution to N(e -C T c ) is from 
the coherence peaks located at momenta p e such that 
e = A y(p £ ), close to the nodal directions, see Fig. lb). 
Vortex scattering is weak for the quasiparticles travel- 
ing nearly parallel to the field, and strong for those with 
the momenta at moderate to large angles with respect to 
H. Consequently, when the field is applied along a nodal 
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Fig. 2. Comparison of the phase diagram for the specific heat as a 
function of Zeeman splitting. In the shaded (unshaded) area C(T, H) 
has a minimum for the field along the nodal (antinodal) direction. 
Here So = (ch/2e) /27t£q with £o = Hvp /2nT c is a measure of the 
orbital critical field. Parameter Z = fiBo/2irT c characterizes the 
strength of the Zeeman coupling. The arrows in the right panel 
indicate onset of the first order transition occurring for Z > 0.35. 

direction, it largely preserves the spectral weight in the 
coherence peaks aligned with the field, while exciting un- 
paired quasiparticles at the other two nodes. In contrast, 
the field applied along the gap maxima induced more low 
energy states, but smears out the coherence peaks at mod- 
erate e, transferring their spectral weight to other energies, 
see Fig. lc). As a result, the DOS anisotropy changes sign 
at a finite e* , which was determined numerically. 

While the results were very suggestive of the d x 2_ y 2 sym- 
metry in CeCoIns , specific details pertinent for this mate- 
rial left some questions open. First, the upper critical field 
in CeCoIns is Pauli-limited [10], and spin splitting of the 
Fermi surface may be important. We considered its effect 
by including the Zeeman term, v z = guBO~ ■ H with g = 2 
into the semiclassical equations, and writing explicitly the 
equations for the two spin components, which are indepen- 
dent. On the other hand, the self-consistency involves spin 
summation. We found that the inclusion of Zeeman split- 
ting modifies the location of the inversion lines in the phase 
diagram, but does not affect the main conclusion, namely, 
the inversion of the anisotropy at moderately low H and 
T, see Fig. 2. The effect of the paramagnetic limiting on 
the phase diagram is stronger at high fields and low tem- 
peratures due to "isotropization" of the upper critical field: 
the orbital limiting field differs for the nodal and antin- 
odal orientation of the field, while the Pauli limiting field is 
direction-independent. Since in practice the measurements 
of the anisotropy are carried out at moderate fields, our 
conclusions remain substantially unmodified. 

Second, the semiclassical analysis suggests that the lo- 
cal curvature of the Fermi surface around the nodal points 
affects the filling of the near-nodal states due to the field: 
for a fixed position and superflow v s (r), the variation of 
the Fermi velocity vp(p) with direction p near the node 
affects the Doppler shift. Moreover, in borocarbide super- 
conductors some of the features of the anisotropy profile 
were suggested to be due to the nested parts of the Fermi 



surface [11] . Consequently, we explored the influence of the 
Fermi surface geometry on our conclusions. 

We consider a model two-dimensional Fermi surface 
given by p 2 F = p 2 . + p 2 + a 2 (p x + p y ), and chose ap F = 2. 
In this case the directions along the axes and p x = p v are 
inequivalent, and therefore the map of the specific heat 
anisotropy in the T-H plane depends on whether the order 
parameter is of the d x 2_ y 2 or d xy type, with nodes along 
[110] or [100]. We consider both cases, and model the order 
parameter by A(p x ,p y ) = A (p 2 x - p 2 y ) and A(p x ,p y ) = 
Ao(2p x p y ) respectively, with p x ,p y at the FS. Therefore 
the nodes are located either in the "flat" (low curvature) or 
in the "corner" (high curvature) part of the Fermi surface. 

While it is not possible to carry out the fully self- 
consistent calculation in field for a FS with no c-axis energy 
dispersion, we showed in Ref.[9j that solving for Aq(T, H = 
0) and then modeling the vortex state by the lowest Landau 
Level with A(T, H) = A (T, 0)^/1 - H/H c2 gives results 
in semi-quantitative agreement with the self-consistent 
calculation for a corrugated FS. This is also seen from com- 
paring the left panels of Fig. 2 (self-consistent, corrugated 
FS) and Fig. 3 (2D FS). We show the results in Fig. 3. It 
is clear from that the location of the anisotropy inversion 
lines is sensitive to the shape of the Fermi surface. 

If the nodes of the gap are at the regions of high FS cur- 
vature, the "semiclassical" regime of minima of the specific 
heat for the field along the nodes is suppressed. Now even 
for the field along the node the near-nodal quasiparticles 
have the velocity at a substantial angle to the field direc- 
tion, and are efficiently scattered by the vortices; hence 
the anisotropy crossover line is pushed to ultra-low fields, 
where we cannot detect it numerically. Over a large part of 
the phase diagram the minima of the specific heat corre- 
spond to the field along the gap maxima. Similarly to the 
rotationally symmetric Fermi surface, the anisotropy is re- 
versed at low T and high fields, but in this region near the 
transition line the amplitude of the oscillations is small, 
and hence is of little experimental relevance. 

In contrast, flattening of the near-nodal Fermi surface 
slightly extends the "semiclassical" region. In this geome- 
try when the field is applied along a nodal direction, the 
Doppler shift nearly vanishes not just at the nodal point 
where k„ is normal to the superfluid velocity, v s , but over 
the entire near-nodal region. Moreover, the scattering of the 
quasiparticles on vortices in this region is small, since their 
velocity is nearly parallel to the field. As a result, semiclas- 
sical results have a somewhat extended regime of validity. 

In CeCoIns the FS is not rotationally symmetric [12] . 
and likely locations of the nodes pass through regions of 
high curvature. Therefore the inversion line is likely to be 
at lower T, iJ, than a cylindrical FS model predicts. 

4. Summary 

We reviewed here recent progress in microscopic calcu- 
lations of the electronic specific heat of unconventional su- 




Fig. 3. Effect of the Fermi surface shape on the anisotropy map of 
the specific heat. Left panel: cylindrical Fermi surface; middle panel: 
nodes in the region of high curvature; right panel: nodes in the 
region of low curvature. In the shaded (unshaded) area C(T, H) has 
a minimum for the field along the nodal (antinodal) direction. 

perconductors under a rotated magnetic field, and its con- 
nection with the experimental efforts to use the method for 
determination of the nodal directions. The approach uses 
the Brandt-Pesch-Tewordt approximation, which is nearly 
exact at high fields, and likely overestimated the effect of 
vortex scattering at low fields. The major difference with 
the results of the semiclassical approximation is the inver- 
sion of the anisotropy of the density of states and the spe- 
cific heat due to the competition between the Doppler en- 
ergy shift and spectral weight redistribution resulting from 
scattering on the vortices. The low-to-moderate T and H 
location of the inversion line makes it of experimental rele- 
vance, and we showed that the field and temperature range 
of this inversion is weakly dependent of paramagnetic ef- 
fects, but is sensitive to the shape of the Fermi surface. A 
more complete analysis of the Fermi surface effects and Zee- 
man splitting, including the behavior of the thermal con- 
ductivity, will be published elsewhere. 

We benefited from enlightening discussions with Y. Mat- 
suda and T. Sakakibara. This research was supported in 
part by the Louisiana Board of Regents. 
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